Evolution of level density step structures from ^^-^^Fe to ^^'^^Mo 
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Level densities have been extracted from primary 7 spectra for ^^'^''Fe and ^®'^^Mo nuclei using 
(^He,a7) and (^He,^He'7) reactions on ^'^Fe and ^^Mo targets. The level density curves reveal step 
structures above the pairing gap due to the breaking of nucleon Cooper pairs. The location of the 
step structures in energy and their shapes arise from the interplay between single-particle energies 
and seniority-conserving and seniority-non-conserving interactions. 
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The group at the Oslo Cyclotron Laboratory (OCL) 
has recently developed a new method (the so-called 
Oslo method) to extract the level density and radia- 
tive strength function from primary 7 spectra [1]. The 
method can be characterized as a further development of 
the sequential extraction method [2,3]. The Oslo method 
has been extensively tested in the rare-earth mass region 
[4-7] and has been successfully extended to the very light 
27,28gj nuclei [8]. The robustness of the Oslo method 
in the mid sd shell vs^here the level density is rather low 
and the single-particle levels are widely spaced prompted 
the present investigations in mass regions around the 
/7/2 shell (the iron nuclei) and around the Jpg shell 
(the molybdenum nuclei). The ^^Fe nucleus was espe- 
cially chosen due to astrophysical interest in this isotope. 
While Ref. [9] discusses the e-process and the direct pro- 
duction of ^^Fe, a modern view in terms of a hierarchy of 
restricted nuclear statistical equilibria (NSE) on different 
time scales predicts the direct production of ^^Fe only 
on the highest level of NSE, namely in the presence of 
weak equilibrium [10]. These rather extreme conditions 
are envisioned only in the core of massive stars (where 
ejection into space is unlikely) and in the core of Type 
la supernovae (where ejection into space seems possible). 
In the equations of NSE, the nuclear partition function, 
i.e., the Laplace transform of the nuclear level density, 
plays a prominent role. The nucleus ^^Mo is of special 
interest in the investigation of the \N — Z\ dependence 
of level densities, since the mass 96 isobars are the only 
ones in the nuclear chart where one can find three differ- 
ent stable nuclei with \N — Z\ varying by eight units from 
^^Zr to ^^Ru. Also the — Z\ dependence of level den- 
sities can have significant astrophysical importance [11], 



since many reactions take place on unstable isotopes and 
therefore most of their reaction cross sections must be 
estimated by Hauser-Feshbach type calculations [12]. 

The purpose of this Rapid Communication is to report 
on experimental level densities of ^^-^'^Fe and ^^'^^Mo nu- 
clei and to provide a schematic explanation of the ob- 
served step structures in these curves. 

The experiments were performed at the MC-35 cy- 
clotron at the OCL using a ~ 2 nA beam of 45 MeV 
■^He particles. The self-supporting targets were isotopi- 
cally enriched to 94.7% and 94.2% and had thicknesses 
of 3.4 mg/cm^ and 2.1 mg/cm^ for ^^Fe and ^^Mo, re- 
spectively. Each experiment ran for approximately five 
days each, and about 200,000 relevant particle-7 coinci- 
dences were recorded in each analyzed reaction channel. 
The charged ejectiles were identified and their energies 
measured in a ring of eight coUimated Si AE-E tele- 
scopes placed at 45° with respect to the beam direction. 
The thicknesses of the front and end detectors were 140 
and 3000 /xm, respectively and shielding against S elec- 
trons was achieved with a 19-/xm-thick Al foil. The dis- 
tance from the target was 5 cm, giving a total solid angle 
coverage of 0.3% of An, and the energy resolution was 
~0.3 MeV over the entire spectrum. The 7-rays were 
detected in 28 collimated 5"x5" Nal(Tl) detectors called 
CACTUS [13] surrounding the target and particle de- 
tectors. The total efficiency of CACTUS was -15% of 
4tt and the resolution was ~6% of the deposited energy 
at 1.3 MeV. In addition, one 60% Ge(IIP) detector was 
used in the setup to monitor the selectivity and populated 
spin distribution of the reactions. Raw a-particle data 
are shown in Fig. 1 for the case of the ^^Mo(^IIe,a)^^Mo 
reaction, where discrete transfer peaks are observed up 
to ~6 MeV in excitation energy indicating the opening 
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of the neutron qqj^ shell. 

From the known Q-value and the reaction kinematics, 
the ejectile energy can be transformed into initial exci- 
tation energy of the residual nuclei. Using the particle-7 
coincidence technique, each 7 ray can be assigned to a 
cascade depopulating a certain initial excitation energy 
in the residual iiuc;lcus. The data arc therefore sorted into 
total 7-ray spectra originating from different initial exci- 
tation energy bins. Every spectrum is then unfolded us- 
ing a Compton-subtraction method which preserves the 
fluctuations in the original spectra and does not intro- 
duce further, spurious fluctuations [15]. Prom the un- 
folded spectra, a primary 7 matrix is constructed using 
the subtraction method of Ref. [16]. The basic assump- 
tion behind this method is that the 7-ray decay pattern 
from any excitation energy bin is independent of whether 
states in this bin are populated directly via the (^He,Q!) 
or ('''He, '"'He') reactions or indirectly via 7 decay from 
higher excited levels following the initial nuclear reaction. 
This assumption is trivially fulfilled if one populates the 
same levels with the same weights within any excitation 
energy bin, since the decay branchings are properties of 
the levels and do not depend on the population mecha- 
nisms. On the other hand, if one populates, e.g., vastly 
different spin distributions within one excitation energy 
bin via the two different population mechanisms, the 7- 
decay patterns should be different as well. As an exam- 
ple of the data analysis discussed in this paragraph, the 
raw, unfolded and primary 7 spectra are shown for the 
57Fe(3He,a7)5'^Fe and "Fe(3He,3He'7)"Fe reactions in 
Fig. 2. 

Finally, the primary 7 matrix is factorized using the 
generalized Brink- Axel hypothesis [17,18]. The original 
hypothesis states that the giant dipole resonance (GDR) 
can be built on every excited state, and that the prop- 
erties of the GDR do not depend on the temperature of 
the nuclear state on which it is built. This hypothesis 
can be generalized to include not only the GDR but any 
type of nuclear excitation and results in the assumption 
that primary 7 spectra originating from the excitation 
energy E can be factorized into a 7-ray transmission co- 
efficient T{E^) which depends only on the 7-transition 
energy Ej and into the level density p{E — Ej) at the 
final energy. This factorization is determined by a least 
fit to the primary 7 matrix, using no a priori assiimp- 
tions about the functional form of either the level density 
or the 7-ray transmission coefficient [1]. An example to 
illustrate the quality of the fit is shown in Fig. 3, where 
for the ^^Mo(^He,^He'7)^^Mo reaction we compare the 
experimental primary 7 spectra from two different initial 
excitation energies to the least fit. Unfortunately, the 
mathematical structure of the relevant equations in the 
least fit does not allow us to find a unique solution 
for the level density and 7-ray transmission coefficient. 
However, it has been shown that all solutions with the 
same x^ can be obtained by the transformation of one 
randomly chosen solution according to [1] 



p{E - Ey) = Acxp[a{E - Ej)]p{E - Ej) (1) 
f{E^) = Bexp{aE^)T{E^). (2) 

The three free parameters A, B, and a have to be de- 
termined to give the physically most relevant solution to 
the least fit using independent experimental informa- 
tion. The most common way is to count the number 
of discrete levels at low excitation energies and to use 
the neutron resonance spacing at _B„ to find values for A 
and a. The remaining parameter B is then determined 
using the average total radiative width of neutron reso- 
nances [6]. Unfortimately, in the case of "''^Fe, there are 
no data on neutron resonances and thus, the informa- 
tion about the level density around i?„ in ^^Fe has to 
be obtained by different means. In order to do so, we 
calculate the level density at i?„ in ^''Fe using a back- 
shifted Fermi-gas expression with the parameterization 
of von Egidy et al. [19], where we apply an additional 
overall re- normalization factor to match the level density 
determined from neutron resonance spacings. Then, we 
use the same level-density formula (including the same 
re-normalization factor but with the other parameters 
for ^^Fe) to calculate the level density at B,, in ^^Fe. 
Using this data point instead of the unknown neutron 
resonance spacing, we proceed in the same way as for the 
other three nuclei [1,20]. 

In Fig. 4 we show the physically most relevant solutions 
for the level densities in 56,57pg g^j^j ^^''^'''Mo from the least 
X^ fit. The most striking feature in those curves are the 
steps starting at 2.9 MeV and 1.8 MeV in ^^Fe and ^'^Fe. 
There are also less pronounced but still statistically sig- 
nificant step structures at 2.0 MeV and 1.2 MeV in ^^Mo 
and ^^Mo. It has been established in the rare-earth re- 
gion that such low-lying step structures are connected 
to the breaking of the first nucleon Cooper pair [21]. It 
is therefore natural to assume that the same is true for 
the lighter mass regions. The additional step at 6 MeV 
in the level density curve of ^^Mo might be due to the 
opening of the 59 /2 shell which is indicated to happen at 
this excitation energy (see Fig. 1). A future systematic 
investigation of Mo nuclei might shed more light on this 
issue. 

The difference in binding energy between the even 
and odd systems is a measure for pairing correlation 
and can be calculated from the three-mass indicator of 
Dobaczewski et al. [22]. This indicator yields 1.2 MeV 
and 0.9 MeV for the Fe and Mo nuclei, respectively, which 
agrees very well with the differences in excitation energy 
for the first steps in the level density curves. Thus we 
know that the steps in neighboring nuclei appear at the 
same 'effective' excitation energy, i.e., the excitation en- 
ergy corrected for the contribution to the binding energy 
due to neutron-pairing correlations. Further, the average 
proton pairing energies are 0.8 MeV and 1.1 MeV for the 
^^Fe and ^^Mo nuclei, respectively. These energies should 
now correspond to the excitation energies of the steps in 
the two odd nuclei. However, the steps are delayed in 
excitation energy by 1.0 MeV and 0.1 MeV for these two 
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nuclei. This might be explained by the fact that one not 
only has to invest the energy to break a proton Cooper 
pair but also at least one of the unpaired protons has to 
be promoted to the next unoccupied single-particle level. 
The spacing of those levels can also be calculated using 
the Dobaczewski throe-mass indicator, giving 2.0 MeV 
and 0.9 McV in the two cases. Clearly, the higher single 
particle spacing for the ^^Fe nucleus leads to a larger de- 
lay in excitation energy for the step structure to appear 
compared to the ^^Mo nucleus, but obviously the exact 
excitation energy for the steps in the level density curves 
cannot be estimated from binding energies, since it will 
depend on the exact location of the Fermi energy within 
the single-particle level scheme. 

Another complication might arise from the effect of se- 
niority non-conserving interactions, which mix configura- 
tions of different seniority and smooth out the step struc- 
tures in the level densities. To investigate this, we have 
performed a model calculation. In the model we assume 
a system of eight particles scattered into an equidistant 
single-particle level scheme with eight doubly-degenerate 
levels. As residual interactions we consider a pairing 
interaction and a seniority non-conserving interaction. 
Thus, the model Hamiltonian is written as 

8 1 * 

^ = ia{a^ - -G ^ ajalajaj 

i=l iij=l 

1 ^ 

--^K ^ Wi.jkia\a]akau (3) 

where and a are Fermion creation and annihilation 
operators and the labels with bars stand for time re- 
versed orbits. The single particle level spacing e, the 
strength G of the pairing interaction and the strength 
K of the seniority non-conserving interaction W are the 
only macroscopic parameters of the model. This model 
with good seniority, i.e., the case of k = has already 
been diagonalized in Ref. [23], and the dotted line in 
Fig. 5 gives the distribution of eigenvalues with excita- 
tion energy, i.e., the exact level density of the model. The 
individual bumps contain mainly levels with the same se- 
niority, thus the step structures in the level density can be 
explained by the consecutive breaking of micleon Cooper 
pairs. However, the experimental data show that in gen- 
eral the steps are much smoother than in such a simple 
model calculation indicating the need for seniority non- 
conserving terms in the Hamiltonian. 

An important contribution to seniority breaking comes 
from quadrupole collectivity, which can change rapidly 
with mass number for the nuclei under study. Strong, 
residual quadrupole-quadrupolc interactions can lead to 
deformation of the nucleus and to a modification of the 
single-particle level spacing. Unfortunately, our sim- 
ple model cannot incorporate a realistic quadrupolc- 
quadrupole interaction; the smoothing of the level den- 
sity, however, should not depend on the details of the 



seniority non-conserving residual interaction. We choose 
therefore to model the quadrupole-quadrupolc interac- 
tion by a random two-body interaction [24] of roughly 
equivalent strength where all the pairing-like terms have 
been set to zero i.e., the Wijki in Eq. (3) are Gaussian 
random numbers of mean zero and width equal to one 
except for the cases oi j = i and k = 1 where W^^ii = 0. 

Since in this more general case Eq. (3) does not have 
good seniority, exact diagonalization of the model Hamil- 
tonian has to be performed within the full model space 
which is computationally very demanding, as one re- 
quires all levels, not just the lowest. For eight particles 
in eight doubly-degenerate states there are 12870 states 
(counting all magnetic substates). The number of lev- 
els (not counting magnetic substates) equals the num- 
ber of states with spin projection = 0, thus yield- 
ing 4900 levels. For the present calculations, we used 
e = 0.25 MeV, G = 0.5 MeV, and k « MeV (pure pair- 
ing) and K = 0.14 MeV (pairing+random interaction). 
In Fig. 5, we show the resulting level densities of the two 
calculations. The resolved bumps with definite seniority 
in the pure pairing case are smeared out by the random 
interaction. The gaps between the bumps are rapidly be- 
ing filled, while the bumps themselves are degraded into 
a step structure quite similar to the step structure seen in 
the experimental level density curves of the iron nuclei. 
The occurrence of the step structure in the calculation is 
actually quite sensitive to the exact choice of the strength 
of the random interaction. A weak random interaction 
(k < 0.1 MeV) does not fill the gaps between the bumps, 
a strong random interaction [k > 0.2 McV) produces 
a more smeared out step structure in the level density 
curve similar to the Mo data, indicating the presence of 
relatively stronger seniority non-conserving interactions 
in the Mo nuclei. The range of values which produces 
good qualitative agreement with the Fe data is about 
0.13 MeV< K < 0.16 MeV. However, due to the simplic- 
ity of the model, no attempt has been made to achieve 
quantitative agreement between the model and the ex- 
perimental data. 

In conclusion, we have presented new experimental 
data on level densities below 5„ in ^^■■'''''Fc and ''*''^^Mo. 
Step structures in the level density curves have been re- 
lated to the breaking of nucleon Cooper pairs. The occur- 
rence of the step structures in energy depends on the ex- 
act location of the Fermi energy in the single particle level 
scheme. The smoothness of the step structures depends 
on the strength of the seniority non-conserving interac- 
tions, e.g., the quadrupole-quadrupole interaction com- 
pared to the pairing interaction and the average single- 
particle level spacing. The more pronounced step struc- 
tures in the level density data of Fe nuclei compared to 
Mo indicate relatively stronger seniority non-conserving 
interactions in the latter case. Our schematic calcu- 
lations where random two-body interactions were used 
to model residual seniority non-conserving interactions 
show that level density curves are very sensitive to the 
relative strengths of the different terms in the model 
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Hamiltonian. The present experimental data therefore 
break new ground for the possible theoretical investiga- 
tion of two of the most important residual interactions 
in atomic nuclei, namely the pairing and the quadrupole- 
quadrupole interaction 
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FIG. 1. A singles a spectra is shown in the upper panel and an a-7 coincidence spectrum is shown in the lower panel for 
the ^^Mo('^He,a)^''Mo reaction. At the neutron binding energy Bn, the coincidence spectrum shows a dip reflecting the lower 
7 multiplicity in the decay from low-lying states in ®^Mo which are populated by neutron emission. The transfer peaks to the 
first three states are well known from the {p,d) reaction [14]. The strong transfer peak at 5.3 MeV is a new discovery and 
indicates the opening of the neutron g9/2 shell at this energy. 
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FIG. 2. Raw, unfolded, and primary 7 spectra from the ''^Fe(^He,Q7)^''Fe reaction at 5 MeV excitation energy (upper panels) 
and from the ®'^Fe(^He,^He'7)^^Fe reaction at 6.2 MeV excitation energy (lower panels). 
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FIG. 3. Experimental primary 7 spectra (data points with error bars) at two different initial excitation energies (indicated in 
the figure) compared to the least fit (solid lines) for the ^^Mo('^He,^He'7)®^Mo reaction. The fit is performed simultaneously 
to the whole primary 7 matrix of which the two displayed spectra are only a small fraction. The first generation spectra are 
normalized to one at each excitation energy. 
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FIG. 4. Experimental level densities for the four nuclei under study. Where available, level density data from neutron 
resonance spacirigs have been added (open triangles). The solid lines are the re-normalized level density parameterizations 
according to von Egidy ct al. [19]. Stop structures in the level densities are marked by arrows. The insert shows the step 
structure for ^''Mo on a linear scale. In the level density of ^®Fe, the bump and the plateau at 0.8 MeV and 2.0 MeV, 
respectively, are due to the first and second excited states. 
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FIG. 5. Model calculation using the Hamiltonian of Eq. (3) with e = 0.25 MeV, G = 0.5 MeV, and ac = MeV (dotted 
line). The level density is given as function of excitation energy. Adding a random two-body interaction with the strength 
K = 0.14 MeV (solid line) results in a step structure quite similar to what is seen in the experimental level density curve of the 
iron nuclei. 
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